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ABSTRACT 

We analyze the quark mass and flavour dependence of the QCD phase transition 
temperature. When the lightest pseudo-scalar meson mass (mpg) is larger than 
2 GeV the critical temperature is controlled by the gluonic sector of QCD alone. For 
smaller values of the lightest meson mass the pseudo-critical temperature decreases 
slowly with mpg. For a large regime of meson masses the pseudo-critical temperature 
of 2-flavour QCD is about 10% larger than in the 3-flavour case. On lattices with 
temporal extent N T = 4 an extrapolation to the chiral limit yields T c = (173±8) MeV 
and (154±8) MeV for 2 and 3-flavour QCD, respectively. We also analyze dynamical 
quark mass effects on the screening of the heavy quark potential. A detailed analysis 
of the heavy quark free energy in 3-flavour QCD shows that close to T c screening 
effects are approximately quark mass independent already for pseudo-scalar meson 
masses mpg < 800 MeV and screening sets in at distances r ~ 0.3 fm. 
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1 Introduction 



One of the most prominent non-perturbative features of QCD is the existence of 
a chiral symmetry restoring and deconfining phase transition at high temperature. 
Details of this phase transition, the quantitative value of the transition temperature 
as well as the order of the transition, will depend on basic parameters of QCD, i.e. 
the number of flavours (up) and the values of the quark masses (m q ). The influence 
of global symmetries of the QCD Lagrangian on the order of the phase transition 
is well understood and predictions in the chiral |l| as well as quenched || limit 
have to a large extent been verified by lattice calculations. It has been shown that 
the chiral transition is first order in 3-flavour QCD 0| and likely to be continuous 
in the 2-flavour case ||. The explicit flavour as well as quark mass dependence of 
the transition temperature itself, however, is not related to universal properties and 
therefore is less well understood. 

Of great interest is, of course, a quantitative determination of the QCD transi- 
tion temperature^ in the limit of light quark masses. The transition temperatures in 
2 and 3-flavour QCD will provide upper and lower limits for the transition tempera- 
ture of QCD with two light and a heavier strange quark mass. Previous calculations 
of the transition temperature in QCD with dynamical quarks have shown that it is 
significantly smaller than the critical temperature determined for the deconfinement 
transition which occurs in the limit of infinitely heavy quarks, i.e. in the SU(3) 
gauge theory ||. The quark mass dependence of the QCD transition temperature, 
however, is not well understood. In particular, we would like to understand in more 
detail what sets the scale for the QCD transition temperature. It is obvious that 
for large quark masses the valence quark dominated hadronic sector will completely 
decouple and thermodynamics will be controlled by the gluonic sector alone. In the 
chiral limit, on the other hand, the light pseudo-scalar mesons are expected to play 
a dominant role. For instance, studies within the context of a— models suggest a 
strong dependence of the transition temperature on the pseudo-scalar meson mass 
0. For realistic quark mass values the transition temperature itself will be com- 
patible with the value of the pseudo-scalar mass and one may expect that the rich 
hadronic resonance spectrum will play an important role in this regime . The im- 
portance of a rapidly (exponentially) rising resonance mass spectrum for the location 
of the transition temperature does become apparent already in the SU(N) gauge 
theories where the transition temperature is rather low (T c ~ 270 MeV) despite the 
fact that the lightest excitations in the low temperature phase are glueballs with a 
mass of about 1.5 GeV ||. In fact, the value of the phase transition temperature of 
SU(N) gauge theories can be well understood in terms of the exponentially rising 
mass spectrum obtained in string models ||. In order to better understand QCD 

a In the following we often will talk about the QCD transition without distinguishing between 
a rapid crossover and a true phase transition. This should be clear from the context. 
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thermodynamics and in particular the critical behaviour close to the transition tem- 
perature a detailed analysis of the quark mass dependence of the QCD transition 
is needed. This problem will be addressed here. In order to be able to discuss the 
quark mass dependence of the QCD transition temperature we also have to address 
the question in how far we can compare transition temperatures when the quark 
masses and/or the number of flavours are varied. We need an appropriate observ- 
able to set the scale for the transition temperature which itself is little influenced 
by these external parameters. In order to reach quantitative conclusions on these 
questions it will be mandatory to get the intrinsic systematic lattice effects, cut-off 
dependence and finite volume effects, under control. We therefore will use here im- 
proved gauge and fermion actions for our analysis of the quark mass and flavour 
dependence of the transition temperature in QCD. 

We also discuss the influence of dynamical quarks on the heavy quark potential, 
or more precisely the heavy quark free energy extracted at finite temperature from 
Polyakov loop correlation functions. In the case of 3-flavour QCD we analyze in 
detail the quark mass and temperature dependence of the screening of the heavy 
quark free energy, which becomes significant already below T c . 

This paper is organized as follows. In the next section we will introduce the 
improved gauge and staggered fermion actions used for our simulations and give 
some details on the simulation parameters. In Section 3 we present our results on 
the flavour and quark mass dependence of the transition temperature and give an 
estimate for the critical temperatures of the chiral transition in massless 2 and 3- 
flavour QCD. In Section 4 we discuss in more detail the quark mass and temperature 
dependence of the heavy quark free energy. We finally give our conclusions in Section 
5. 

2 Thermodynamics with an improved staggered 
fermion action 

The use of improved gauge and fermion actions is mandatory in thermodynamic 
calculations which generally are performed at quite large values of the lattice spac- 
ing, a, because one has to use quite small values for the temporal extent of the 
lattice. The finite cut-off effects become particularly striking in the calculation of 
bulk thermodynamic observables. In fact, it has been shown just in these cases that 
improved actions can be particularly efficient and reduce the cut-off dependence of 
thermodynamic observables significantly ||. 

In the pure gauge sector the tree level, Symanzik improved action yields a sat- 
isfactory description of the high temperature ideal gas limit of QCD already on 
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lattices with temporal extent N T = 4 P0[ . Similarly good results can be achieved 
in the fermion sector with improved staggered fermions where three-link terms are 
added to the ordinary one link discretization of the first order derivatives appearing 



in the fermion Lagrangian [11, 12 



In the present analysis we will use such an improved action, which in addition 
to the standard Wilson plaquette term and the standard 1-link staggered fermion 
action also includes the planar 6-link Wilson loops and bended 3-link terms which 
improve the rotational symmetry of the staggered fermion action [ 12 1 . The partition 
function for QCD with n/ quark flavours reads 



Z(T, V) = [ ]JdU,(x)e-^ J] ( f Y[dx x d Xx e 

J x,n q=l \ J x 



-S F {m q ) 



1/4 



(2.1) 



where x labels the sites on a hypercubic lattice of size N% x N T . The gauge (Sq) 
and fermion (Sp) actions are given by 
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Here tj^x) = (— l) a: o+--+ a: M-a denotes the staggered fermion phase factors. Further- 
more, we have made explicit the dependence of the fermion action on different quark 
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flavours q, and the corresponding bare quark masses m q , and give an intuitive graph- 
ical representation of the action. The tree level coefficients cf and of appearing in 
Sf have been fixed by demanding rotational invariance of the free quark propagator 
at 0(p 4 ) ( "p4-action" ) Jll|. In addition the 1-link term of the fermion action has 
been modified by introducing "fat" links |H| with a weight u = 0.2. The use of fat 
links does lead to a reduction of the flavour symmetry breaking close to T c and at 
the same time does not modify the good features of the p4-action at high temper- 
ature, i.e. it has little influence on the cut-off dependence of bulk thermodynamic 
observables in the high temperature phase 111], |i~2|. Further details on the definition 



of the action are given in [11 



For simplicity, in the following we will refer to the p4-action with a fat 1-link 
term combined with the tree level improved gauge action as the p4-action. 

Numerical simulations with dynamical fermions have been performed using the 



Hybrid R algorithm [I4| . This algorithm introduces systematic errors which depend 
quadratically on the step size 5t used for the integration of molecular dynamics 
trajectories. We have verified on a small 8 3 x 4 lattice that the systematic errors 
with our improved action are similar in magnitude to those found for the standard 
action ||15|| . In our simulations we use a step size St = min{0.4 m q , 0.1}. The length 



of the molecular dynamics trajectories has been chosen to be r = 0.8. 



3 Critical Temperature in 2, 2+1 and 3 Flavour 
QCD 

In this section we want to discuss our results for the pseudo-critical temperatures 
obtained for 2 and 3-flavour QCD with different quark masses. In addition to the 
results obtained with identical quark mass values for all flavours we also will present 
the result of a calculation for (2+l)-flavour QCD, i.e. for QCD with two light 
quark flavours (m U]£ z = 0.1) and one heavier flavour (m s = 0.25). Of course, these 
bare quark mass parameters receive a multiplicative renormalization. For an easier 
comparison with the better known parameters of the unimproved staggered fermion 
action we note that the bare quark parameter, m q , of the p4-action corresponds to 
about twice the bare quark mass in the unimproved staggered action. 

The calculations with degenerate quark masses cover a range of pseudo-scalar 
meson mass from about 250 MeV up to the quenched limit of infinitely heavy mesons. 
We will analyze the quark mass dependence of the pseudo-critical temperature and 
discuss its extrapolation to the chiral limit. This involves two steps: First of all one 
has to determine the pseudo-critical couplings, /3 c (m q ), on finite temperature lattices 
of size x N T . One then has to calculate physical observables at zero temperature 
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n f = 2 


m u ,d 


iV 3 x N T 


#/? 


#iter. 




0.025 


8 3 x 4 


3 


3200 


3.542 (7) 


0.05 


8 3 x 4 


4 


3700 


3.585 (11) 


0.10 


16 3 x 4 


10 


36300 


3.646 (4) 


0.20 


8 3 x 4 


6 


6300 


3.778 (12) 


Tlf = 2 + 1 


m u4 /m s 


iV 3 x N T 


#/? 


#iter. 


ft 


0.10/0.25 


16 3 x 4 


10 


31550 


3.543 (2) 



Table 1: Simulation parameters for the determination of the critical temperature 
for 2 and 2 + 1 flavours. The largest lattice size simulated, the number of (3 values 
and iterations as well as the finally determined pseudo-critical couplings are given. 

which allow to extract the lattice spacing at (3 c {m q ) through which the temperature 
is fixed, T c = l/N T a(/3 c ). 



3.1 Critical couplings on 16 3 x 4 latt ices 

The finite temperature calculations have been performed on lattices with temporal 
extent N T = 4 and spatial lattice sizes N a = 8, 12 and 16. For each value of the 
quark mass we have performed simulations at several values of the gauge coupling 
p. Details on the simulation parameters and the statistics collected on our largest 
lattices are given in Table [I]. For the determination of the pseudo-critical couplings 
we calculated the average Polyakov loop, 
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(3.1) 



and used a noisy estimator for the chiral condensate which is based on 25 random 
vectors FL, 
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As usual we have introduced here the fermion matrix M, which is defined through 
the matrix form of Eq. p.3| , i.e. Sp = Y^q=i X M(m q ) x- At each value of the quark 
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mass the data from runs at all /3-values have been used in a Ferrenberg-Swendsen 
reweighting [TO to determine the maxima in the related Polyakov-loop and chiral 
susceptibilities^], 
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(3.3) 



as well as the susceptibility of the gauge action, 

Xs = N*N T ((S G ) - (S G ) 2 ) . (3.4) 

At intermediate values of the quark mass the transition to the high temperature 
phase of QCD is not related to any form of critical behaviour in thermodynamic 
observables. Nonetheless, the maxima of the susceptibilities define pseudo-critical 
couplings, /3 c (m q ), for the transition to the high temperature phase of QCD, which 
manifests itself in a sudden change of thermodynamic observables in a small tem- 
perature interval^]. For small (vanishing) and large values of the quark mass the 
transition, however, becomes a true phase transition in the infinite volume limit. In 
these cases the pseudo-critical couplings give finite volume estimates for the location 
of the critical point in the thermodynamic limit. 

Previous analyses of the volume dependence of the pseudo-critical couplings per- 
formed with unimproved gauge and staggered fermion actions j|, [15] have shown 
that finite size effects are small at intermediate values of the quark mass. Moreover, 
the pseudo-critical couplings seem to show little dependence on the observable used 
to define them. This is also supported by our analysis of finite volume effects in 
simulations with the p4-action. In the case of 3-flavour QCD we have determined 
the pseudo-critical couplings /3 c (m 5 ) for 10 different values of the bare quark mass 
in the interval 0.01 < m q < 1.0 on lattices of spatial extent N a = 8, 12 and 16. The 
results of this analysis are summarized in Table ^| Within the statistical errors of 
the present analysis we do not observe any clear volume dependence at intermediate 
values of the quark mass and even not for the smallest value, m q = 0.01, considered 
here. This, however, is no longer true for large masses, where the onset of the first 
order deconfining transition leads to large correlation lengths. This is reflected in 
a strong volume dependence of the peak height of the susceptibilities and a visible 
volume dependence of the pseudo-critical couplings. Note that the location of the 

b Note that we only calculate the disconnected part of the complete chiral susceptibility, Xz^t, ~ 
8 2 \nZ/d(m q ) 2 . 

c It recently has been suggested that the sudden change in behaviour at a certain pseudo-critical 
temperature may be characterized in terms of a non-thermal, percolation phase transition [|l7| . 
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pseudo-critical couplings (3 C 
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N a 


#/? 
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max(xs) 


max(xL) 


max (Xiptp,disc) 
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10300 


3.2865 (98) 


3.2887 (65) 


3.2841 (57) 
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4440 


3.2784 (89) 


3.2804 (37) 


3.2806 (20) 


025 


8 


8 


14000 


3.3265 (93) 


3.3318 (24) 


3.3281 (25) 


16 


6 


8200 


- 


3.3428 (159) 


3.3194 (21) 


05 


8 


11 


16250 


- 


3.4018 (35) 


3.3930 (201) 


16 


4 


4500 


3.3866 (47) 


3.3950 (150) 


3.3877 (26) 


0.10 


8 


12 


25100 


3.4637 (34) 


3.4864 (25) 


3.4842 (151) 


16 


9 


15350 


3.4752 (10) 


3.4856 (138) 


3.4756 (14) 


0.20 


8 


16 


39250 


3.5184 (462) 


3.6216 (27) 


3.6204 (32) 


12 


5 


18700 


3.5986 (25) 


3.6020 (32) 


3.6004 (18) 


16 


8 


29650 


3.5998 (14 ) 


3.6017 (17) 


3.6000 (33) 


0.40 


8 


11 


36600 


3.7756 (64) 


3.7848 (23) 


3.7808 (27) 


12 


3 


11700 


3.7810 (24) 


3.7781 (17) 


3.7797 (16) 


16 


5 


16250 


3.7653 (16) 


3.7719 (33) 


3.7746 (40) 


0.60 


8 


12 


33100 


3.7983 (159) 


3.8865 (17) 


3.8806 (12) 


12 


4 


16500 


3.8718 (47) 


3.8770 (16) 


3.8798 (38) 


16 
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41750 


3.8210 (233) 


3.8766 (6) 


3.8764 (6) 


0.70 


16 
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24000 


3.7983 (159) 


3.8865 (17) 


3.8806 (12) 


0.80 


16 


3 


25300 


3.9358 (18) 


3.9377 (10) 


3.9368 (14) 


1.00 


8 


8 


21750 


3.9680 (6) 


3.9686 (27) 


3.9663 (11) 


12 
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24520 


3.9717 (39) 


3.9712 (22) 


3.9711 (15) 


16 
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101300 


3.9771 (15) 


3.9778 (10) 


3.9808 (16) 


oo 
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129000 


4.0554 (27) 


4.0745 (17) 
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12 


14 


154900 


4.0708 (27) 


4.0708 (6) 


- 


16 


8 


157500 


4.0719 (7) 


4.0715 (5) 


- 


24 


3 


75200 


4.0724 (4) 


4.0722 (4) 




32 


1 


46800 


4.0733 (4) 


4.0729 (3) 




oo 






4.0732 (4) 


4.0730 (3) 





Table 2: Pseudo-critical couplings of the p4- action calculated on lattices with 
spatial extent N a = 8, 12 and 16 and temporal extent N T = A for 3-flavour QCD 
with different quark masses. Errors on the maxima of the susceptibilities have been 
determined from a jackknife analysis. 
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peak of the action susceptibility xs shows the strongest volume dependence. How- 
ever, on the largest lattice (16 3 x 4) it also agrees within statistical errors with all 
the other observables used to locate the pseudo-critical points. In particular, we do 
not find any systematic volume dependence in the peak location of xl and Xm disc 
for quark masses m q < 0.6. The pseudo-critical couplings taken from the largest 
lattice which have been used for the determination of the critical temperature are 
also given in Table [p. 

The extrapolation of the pseudo-critical couplings to the chiral limit will depend 
on the order of the phase transition which, in accordance with present knowledge 
from lattice calculations, is expected to be second order for nf = 2 and first order 
for rif > 3 If the transition is first order in the chiral limit it will remain to be 
first order also for small values of the quark mass. The shift in the critical coupling 
will then depend linearly on the quark mass with a slope controlled by the zero 
quark mass discontinuities in the order parameter and gauge action. From a Taylor 
expansion of the free energy below and above the critical point one obtains^, 

Pc(m q ) = &(0) + ^y™, , n f >3 . (3.5) 

In the case of a second order transition, however, the diverging correlation length at 
the critical point will lead to a stronger variation of the critical couplings with the 
quark mass. From the non-analytic structure of the singular part of the free energy 
one finds that the critical couplings will scale like, 

&K)=&(0) + c-K) 1//35 , n f = 2 . (3.6) 

Here 1/(38 is a combination of critical exponents which for the case of 2-flavour QCD 
are expected to be those of the three dimensional 0(4) symmetric spin models, i.e. 
in the continuum limit one expects to find l/[38 ~ 0.55. Of course, the regular part 
of the free energy will also add a linear dependence on m q to this leading order 
behaviour. 

In Figure [I] we show our results for f3 c (m q ) in 2 and 3-flavour QCD obtained for 
small and intermediate values of the quark mass with the p4-action. Furthermore, 
we show results from calculations with the standard Wilson gauge and standard 
staggered fermion action |TS|, [H5], ^U], ^T], Q 

The similarity in the quark mass dependence of f3 c {m q ) seen in Figure [I] for the 
2 and 3-flavour cases suggests that for nf=2 also the subleading corrections, linear 

d This follows the same line of arguments given in Ref. for the quark mass dependence of 
the deconfinement transition in the limit of heavy quarks. 
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Figure 1: Pseudo-critical couplings of 2 and 3 flavour QCD versus bare quark 
mass calculated on lattices of size 16 3 x 4 with the standard Wilson plaquette and 
standard staggered fermion action (std) and the improved gauge and fermion action 
(p4) described in section 2. The data for the standard action have been shifted 
by A/3 = —1.5. Lines show linear fits to the 3-flavour data for m q < 0.05 and 
"0(4)+linear" fits for the 2-flavour data with m q < 0.2. 



in m q are still important in the quark mass regime explored by present simulations. 
We thus have extrapolated the critical couplings to the chiral limit by assuming the 
0(4) symmetric ansatz given in Eq. |3.6| and allowing for a subleading term linear in 
m q . This yields for 



n i 



2: 



A(o) 



5.232 ±0.002 
3.48 ± 0.03 



standard gauge and fermion action 
improved gauge and fermion action 



(3.7) 



The result for the standard staggered action is consistent with fits solely based on 
the ansatz given in Eq. [0] with free exponent 1/(35 [T5|, E3| as well as fits in which 
the exponent has been fixed to its 0(4) value, i.e. 1/(35 = 0.55. 

In the 3-flavour case, on the other hand, the linear dependence on m q is evident 
for small quark masses and linear fits with a small \ 2 / (d.o.f.) can be performed for 
m q < 0.05. From these fits we find for the critical couplings of the 3-flavour theory 
in the chiral limit, 
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n f = 3 



a (q\ — J 5-083 ± 0.004 standard gauge and fermion action , , 

\ 3.258 ± 0.004 improved gauge and fermion action 

The quark mass dependence of (3 c (m q ) is significantly different in the case of 
improved and unimproved actions. From the fits of the critical couplings we find, 

A(^) f 0.49 ± 0.01 standard action 

(3.9) 



A(S G ) ■ N T I 0.70 ± 0.04 improved action 

However, in order to relate these slope parameters to a physical observable, the ratio 
of the discontinuity in the chiral condensate and the latent heat Ae, on still has to 
determine a normalization constant Z m (/3 C ) [^4] that turns (ipip) into a renormaliza- 
tion group invariant observable, 

Al ~ A(S G ) ■ N T Zm {Pc) ■ (3 - 10) 

This may partly account for the observed difference between the standard and im- 
proved actions. 



3.2 Setting the scale at zero temperature 

In order to relate the pseudo-critical couplings to physical temperatures and to 
perform the extrapolation to the chiral limit we have performed zero temperature 
simulations on lattices of size 16 4 . In addition to pseudo-scalar (mpg) and vector 
(my) meson masses we also have calculated the string tension (a) from large Wilson 
loops. 

The spectrum calculations performed at f3 c (m„) follow the conventional approach. 



They are based on wall source operators as described e.g. in Ref. [[25J. Masses have 
then been extracted from the exponential decay of correlation functions constructed 
from wall source and local sink operators. Also the calculation of the string tension 



is based on a standard approach. We use smeared Wilson loops |26] from which we 
extract a potential which then is fitted to an ansatz which includes a Coulomb and 
a linearly rising term. The notion of a string tension in the presence of dynamical 
quarks does, however, require a bit more explanation. In the presence of dynam- 
ical quarks the heavy quark potential is no longer linear at large distances as the 
spontaneous creation of quark anti-quark pairs from the vacuum will lead to string 
breaking. Nonetheless, it has been found that the "potential" extracted from Wilson 
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loops does not show evidence for string breakingF] even at distances of up to about 
2fm [E7|. Wilson loops thus can still be used to extract a "string tension" from 



its large distance behaviour. As will become clear in the following a comparison 
with quenched and partially quenched spectrum calculations shows that this string 
tension shows only little quark mass and flavour dependence. 

The meson masses and the string tension calculated at j3 c (m q ) for 2 and 3-flavour 
QCD are given in Table |3|. Here we also give the result from a simulation of (2+1)- 
flavour QCD. Furthermore, for rif = 3 we have performed calculations with bare 
quark mass m q = 0.1 at various f3- values that cover the region of critical couplings 
obtained with small quark masses. We will use these data in the next section to 
determine the vector meson mass (m p ) in the chiral limit at (3 c (m q = 0). 

In order to quantify the flavour and quark mass dependence of the QCD transi- 
tion temperature it is, of course, important to understand to what extent the zero 
temperature observables used to set the scale are dependent on m q and n/. Clearly 
the hadron masses will depend on the bare quark masses m q . With decreasing quark 
masses one also might expect some n /-dependence resulting from sea quark contri- 
butions to the hadron masses. Quenched spectrum calculations, however, suggest 
that the latter effect is small compared to the explicit valence quark mass depen- 
dence p9 |. The string tension, on the other hand, is expected to show a much weaker 



dependence on the quark masses. This is supported by the fact that in the quenched 
limit, i.e. in the pure SU(3) gauge theory the value for the string tension deduced 



from the ratio (\fa/m p ) quenched hi the limit of vanishing valence quark mass [30 



is consistent with phenomenological values used as input in heavy quark spectrum 



calculations 31 



In Figure |2| we show the ratio ^/a/my as a function of (mps/my) 2 - 111 addition 
to the results obtained from simulations with dynamical quarks at various values 
of the bare quark mass, we also show in this figure results obtained in quenched 
QCD [Q and from a partially quenched analysis. For the latter we have performed 
a spectrum analysis in 3-flavour QCD with valence quarks of mass m q = 0.02 and 
0.05 on gauge field configurations generated with dynamical quarks of mass m q = 0.1 
at the critical coupling, [3 c {m q = 0.1) = 3.475. Results for the partially quenched 
meson masses obtained in this case are given in Table The square root of the 
string tension in units of the vector meson mass in the limit of vanishing valence 
quark mass is obtained as 



v/ct _ f 0.552 ± 0.013 , quenched [30] . y 

m p 1 0.532 ± 0.018 , partially quenched for m q = 0.1, nj = 3 ' 



c String breaking has been found to occur at non-zero temperature where the heavy quark free 
energy can be extracted from Polyakov loop correlation functions (|^] . We will discuss this in more 
detail in Section 4. 
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0.194 (10) 


0.05 


3.395 (15) 


0.303 (13) 


0.706 (6) 


1.320 (35) 


0.454 (10) 


0.189 (6) 


0.10 


3.475 (2) 


0.283 (11) 


0.967 (1) 


1.415 (15) 


0.470 (9) 


0.177 (2) 


0.20 


3.602 (3) 


0.248 (4) 


1.322 (2) 


1.608 (9) 


0.502 (4) 


0.155 (1) 


0.40 


3.772 (4) 


0.189 (4) 


1.814 (4) 


1.985 (10) 


0.575 (6) 


0.126 (1) 


0.60 


3.877 (2) 


0.176 (3) 


2.210 (4) 


2.347 (12) 


0.596 (5) 


0.107 (1) 


1.00 


3.978 (2) 


0.154 (2) 


2.838 (6) 


2.979 (15) 


0.637 (4) 


0.084 (1) 



Table 3: The string tension, pseudo-scalar masses and vector meson calculated at 
the critical couplings which are given in the second column. In the last two columns 
we give the critical temperature in terms of the string tension and vector meson 
mass. 



This clearly is also consistent with an extrapolation of the ratio ^fajmy to the chiral 
limit which yields y/a/m p ~ 0.50. Calculated at non-zero values of the bare quark 
mass this ratio, however, does show a strong dependence on m q . 

We thus find similar results for the ratio y/a / m p in the quenched and in the chiral 
limit. This suggests that the string tension extracted from Wilson loops as well as 
partially quenched meson masses depend only weakly on the number of flavours and 
the value of the dynamical quark mass used in a numerical simulation. They thus 
are suitable observables to set a scale for other physical observables which may be 
sensitive to these external parameters of QCD. 
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m q 


mpga 




0.02 


0.448 (2) 


1.08 (1) 


0.05 


0.695 (2) 


1.23 (2) 



Table 4: Partially quenclied spectrum results in 3-flavour QCD performed at 
Pcim-q = 0.1) on 16 4 lattices. Results for m q = 0.1 and the corresponding value 
of the string tension at the critical coupling are given in Table [| 
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Figure 2: Square root of the string tension in units of the vector meson mass as a 
function of the ratio of pseudo-scalar and vector meson masses. Shown are results for 
2, 2+1 and 3 flavour QCD obtained with the improved gauge and staggered fermion 



action given in Eqs. ^]2| and |2.3| . Also shown are results from calculations with 
unimproved staggered fermions pi| . The black squares show results from quenched 
calculations and partially quenched calculations with a sea quark mass of m q = 0.1. 



3.3 Quark mass dependence of the transition temperature 

In the previous subsections we have discussed the calculations performed to deter- 
mine the pseudo-critical couplings for the QCD transition to the high temperature 
plasma phase and presented zero temperature calculations of meson masses and the 
string tension which are needed to set the scale for the transition temperature. The 
resulting pseudo-critical temperatures expressed in units of the square root of the 
string tension as well as in units of the vector meson mass are given in Table ^ 
and shown in Figure [| We note that T c j \/a and T c /my show a consistent flavour 
dependence. The transition temperature decreases slightly when the number of 
flavours is increased. The effect seems to be somewhat stronger in the latter ratio, 
which reflects the relative flavour dependence of the two observables used to set 
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Figure 3: The transition temperature versus (mps/my) 2 for 2 and 3-flavour QCD 
obtained from calculations with the p4 action on lattices with temporal extent N T = 
4. Also shown are results from calculations using unimproved gauge and staggered 
fermion actions [pOfl . 



the physical scale (see Figure ||). The quark mass dependence or, correspondingly, 
the dependence on the meson mass ratio mpg/my, is, however, completely differ- 
ent for T c /yfa and T c /my. In fact, the decrease of T c /my with increasing quark 
mass mainly reflects the strong quark mass dependence of my itself and does not 
represent the physical quark mass dependence of the transition temperature. The 
ratio T c j yfo ', on the other hand, shows the expected quark mass dependence; with 
decreasing quark mass the transition temperature drops from the pure gauge value, 
(T c /- v /a)g U (3) — 0.635 to a value (T c / '-y/o^cWrai — (0-4 — 0.45). This is in accord with 
the physical picture that with decreasing quark mass hadronic degrees of freedom 
become lighter and thus get more easily excited in a thermal heat bath. They then 
can contribute to the overall energy and particle density of the system and thus can 
trigger the onset of a phase transition already at a lower temperature. Figure |3| thus 
underlines the importance of a physical observable for setting the scale of T c which 
itself is quark mass independent or at least only shows a weak dependence on m q . 
Only in this way we can hope to make contact with phenomenological models that 
discuss the quark mass dependence of T c . We have argued in the previous subsection 
that the string tension and partially quenched meson masses fulfill this requirement. 
We will mainly use here the former observable to set a physically sensible scale for 
T c . Of course, an analysis of T c in units of a hadron masses is still of particular 
importance in the chiral limit as this provides an experimentally well determined 
input whereas y/a is only indirectly accessible through heavy quark phenomenology. 

In Figure [3] we also have shown results for T c j yfa obtained from simulations with 
unimproved gauge and staggered fermion actions on lattices with temporal extent 
iV T = 4. They differ systematically from the results obtained from simulations with 
improved actions, which indicates that we have to consider the influence of finite cut- 



15 



off effects on our results. In the case of unimproved staggered fermions simulations 
of 2-flavour QCD have also been performed on lattices with larger temporal extent, 
i.e. at smaller lattice spacing ||34|| . These calculations indeed lead to larger values 



for T c j yfa which are consistent with our results obtained with improved fermions. It 
thus seems that the simulations with improved gauge and staggered fermion actions 
lead to a sizeable reduction of cut-off effects. 

We have good control over the transition temperatures down to meson mass 
ratios nips/my — 0.4, which corresponds to pion masses of about 370 MeV when 
one uses the string tension to set a physical scaleQ. The results obtained for T c / ^fa in 
2 and 3-flavour QCD indicate that the transition temperature is lower in the latter 
case. From a linear interpolation of our data at intermediate values of the quark 
masses we find 

A(T c /y/a) = f-p) - (-fL) ~ 0.03 for 0.4 < m PS /m v < 0.70 (3.12) 



<J / nf =2 \V cr /n/=3 



Similarly we find A(T c /mv) — 0.015. This corresponds to a temperature difference 
of approximately (10 — 15) MeV. 



3.4 Transition temperature in the chiral limit 

As pointed out in section 3.1 the pseudo-critical couplings will scale differently close 
to the chiral limit of 2 and 3-flavour QCD. This, of course, translates into a differ- 
ent scaling of the critical temperatures as a function of the mass of the Goldstone 
particles, 

T c ( m7T )-T c (0)~\ q 2 71 ' f ~l • (3-13) 
[m q ~ m% , rif > 3 

However, from Figure ^| it is apparent that in the mass regime currently explored by 
our simulations the dependence of T c on the Goldstone mass is quite similar for 2 
and 3-flavour QCD. As also observed in the analysis of the pseudo-critical couplings 
discussed in section 3.1 subleading corrections will be of particular importance for 
the chiral extrapolation in 2-flavour QCD. Subleading contributions linear in m q 
will become relevant and will influence the dependence of T c as well as that of the 
pseudo-scalar mass on the quark masses. In fact, we find that in a large interval 
of pseudo-scalar meson masses T c depends linearly on mps- This is apparent from 



f Whenever we use the string tension to convert to physical units we use \ fa = 425 MeV, which 



is consistent with the quenched and partially quenched results given in Eq. 3.11 and corresponds 
to \falm = 0.55. 
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Figure 4: The transition temperature for 2 and 3 flavour QCD in units of the 
string tension versus mps/v^" obtained with standard (std) and improved (p4) 
staggered fermions on lattices with temporal extent N T = 4. The hatched band to 
the right of the figure denotes the quenched result, the vertical line to the left is the 
physical mps/ ' \fo value. 

Figure [|. The line shown there is obtained from a fit to the 3-flavour data in the 
mass interval, 0.85 < mps/y/cr < 5.3, i.e. for 360 MeV < mps < 2.3 GeV. This 
gave 

(^)K S ) = 0.40 (1) + 0.039 (4) . (3.14) 

In view of the quite similar functional dependence of T c on mps for 2 and 3-flavour 
QCD it seems that this is not indicative for any form of chiral behaviour, although an 
almost linear dependence of T c on mps is just what one would expect to happen in 
the case of 2-flavour QCD {2/(35 = 1.1). We also note the rather weak dependence 
of T c on mps. This suggests that the value of T c is dominantly controlled by a 
large number of resonances at intermediate mass values which is not so sensitive to 
modifications of the mass of the lightest states in the spectrum. 
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3.4.1 Chiral transition in 3-flavour QCD 



In order to get some control over systematic errors resulting from subleading quark 
mass dependences in the various observables used in our analysis we have performed 
extrapolations to the chiral limit in different ways: 



i) We extrapolate T c j 'y/a using linear fits in (m PS /mv) 2 , (jnpsl \f&) 2 an d m q 
for mps/m v < 0.7. This yields T c j yfd = 0.407(15), 0.419(15) and 0.418(16), 
respectively. 

ii) We repeat (i) for T c /m v . This yields T c /m p = 0.204(5), 0.200(3) and 0.200(3). 

iii) We use the vector meson masses calculated at (3 c (m q = 0.05) and f3 c (m q = 
0.025) and extrapolate them to the critical point f3 c {m q — ► 0) = 3.258(4). 
This yields T c /m p = 0.199(4). 



We note that (i) yields results consistent with the fit linear in mps given in Eq. |3. 14j . 
We also stress that the results obtained from (i) and (ii) imply y/a/m p = 0.49(2) 
for the ratio of the two zero temperature observables used to set the scale. This is 
consistent with the results shown in Figure The above extrapolations also are 
consistent with the linear interpolation given in Eq. |3.14| which is shown in Figure |j. 



Combining the results of the various fits performed we estimate for the critical 
temperature of 3-flavour QCD in the chiral limit 

T 

— = 0.20 ±0.01 «-> T c = (154 ± 8) MeV . (3.15) 

m p 

We note that this result has been obtained on lattices with temporal extent N T = 4. 
Although we expect that cut-off effects are small in our calculation with improved 
gauge and fermion actions a systematic analysis on lattices with larger temporal 
extent will be needed to perform a controlled extrapolation to the continuum limit. 
However, we estimate that the systematic error resulting from remaining cut-off 
effects is of the same size as the statistical error quoted above. 



3.4.2 Chiral transition in 2-flavour QCD 



As indicated in Eq. |3.13] , close to the chiral limit we expect to find a universal scaling 
behaviour of T c that is characteristic for the universality class of the second order 
phase transition. For the quark mass or pion mass values accessible in present com- 
puter simulations subleading corrections will, however, still be important. Moreover, 
the explicit quark mass dependence of the zero temperature observable used to set 
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the scale has to be taken into account. This is most drastically seen in the ratio 
T c /my, which decreases with increasing pion mass although one clearly expects that 
the critical temperature in physical units increases with increasing mass of the light- 
est particles in the heat bath. This is correctly reflected in the ratio T c j y/a. In the 
entire quark mass regime analyzed so far the quark mass dependence of the ratio 
T c /my, however, is dominated by the explicit quark mass dependence of my. An 
appropriate ansatz for the quark mass dependence of T c j my is 

(^-) (m g ) = (0) + a mf 35 + b m q , (3.16) 

\myj \m v J y 

where one should find a > and b < 0. Thus, also this ratio will ultimately approach 
the chiral limit from above. With the presently used quark masses, however, the 
subleading term, linear in m q , apparently still gives the dominant contribution. For 
2-flavour QCD the ratio T c /my thus clearly does not yet reflect the expected leading 
order quark mass dependence of T c . 

In performing the extrapolation to the chiral limit we followed the strategy out- 
lined in the previous section. For the case of the ratio T c j yfa we have performed 
extrapolations linear in m q as well as assuming the 0(4) universal behaviour ~ m"' 55 . 
In the case of T c /my, however, we only use the linear extrapolation. From this we 
find 

— = 0.225 ±0.010 , -5= = 0.425 ±0.015 (3.17) 

We note again that this implies for the ratio of the two observables used to set 
the scale, \fojm, p = 0.53(3). Using the rho meson mass to convert the critical 
temperature to physical units we find 

T c = (173 ± 8) MeV . (3.18) 

This result is consistent with earlier estimates obtained from calculations with the 
standard staggered fermion action || and also with a recent calculation performed 
with clover improved Wilson fermions and a renormalization group improved gauge 
action . We stress again that the errors quoted above are statistical only. Similar 
to the nj = 3 case we expect that systematic errors are of similar size. 

Let us finally comment on the critical temperature in the physically most inter- 
esting case of (2±l)-flavour QCD. The simulations performed in this case still use 
too heavy quarks in the light quark sector, i.e. m u ^ = 0.1. However, even in this 
case we find that the pseudo-critical temperature is closer to that of 2-flavour QCD 
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Figure 5: Quark mass dependence of the heavy quark potential for three flavour 
QCD below the transition temperature at a temperature T ~ 0.97 T c . The band 
of lines gives the Cornell-potential in units of the square root of the string tension, 
V{r) 1 1 -Jo = — a/r^Jo+r-Jo with a = 0.25±0.05.The gauge couplings corresponding 
to the different quark mass values are (3 = 4.06, 3.97, 3.86, 3.76, 3.59, 3.46, 3.38, 
3.32 (from top to bottom). 

with the same quark mass than to the 3-flavour case. We thus expect that the crit- 
ical temperature of 2-flavour QCD provides a good approximation to the physically 
realized case of two light and a heavier strange quark. 



4 The heavy quark free energy 

We have discussed in the previous sections the influence of dynamical quarks on the 
QCD transition temperature. In addition to the quark mass dependence of T c we also 
expect that the presence of dynamical quarks of finite mass in the thermodynamic 
heat bath will lead to large qualitative and quantitative modifications of the heavy 
quark potential| V(r,T). In the pure gauge {m q — > oo) limit the heavy quark 
potential is strictly confining for all temperatures below T c . For any finite quark 
mass value, however, string breaking will occur and the potential will approach a 
constant value for r — > oo. 

s We stick here to the commonly used notion of a heavy quark potential, although at non-zero 
temperature we are actually dealing with the heavy quark free energy |35| . 
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Figure 6: Temperature dependence of the heavy quark potential for three flavour 
QCD with a quark mass m q = 0.1. The band of lines gives the Cornell-potential 
in units of the square root of the string tension, V{r)/^/a = —ajr^fo^ r^fo with 
a = 0.25 ± 0.05. The gauge couplings corresponding to the different temperatures 
are (3 = 3.25, 3.30, 3.35, 3.40, 3.43, 3.45, 3.46, 3.50, 3.54. The string tension values 
used to set the scale are based on the interpolation formula given in Ref. 12. The 
potentials have been normalized at short distances such that they agree with the 
zero temperature Cornell potential at r = 1/4T. 

While it seems to be difficult to detect string breaking in calculations of the 
heavy quark potential at zero temperature using Wilson loops, string breaking does 
quite naturally occur at finite temperature, where the heavy quark potential can be 
extracted from Polyakov loop correlation functions 



V gq (r,T) 



— In (L(x)L'(y)) + const. , with \x — y\ 



(4.1) 



In fact, string breaking is directly related to the no n- vanishing of the Polyakov loop 
expectation value at any temperature and finite quark mass, 



lim 

r— >oo 



V qq {r,T) 



T 



-2 In | (L) | + const. 



(4.2) 



First results on string breaking at finite temperature have been reported in 2S 



In order to analyze the quark mass dependence of the heavy quark potential 
and its structure at large distances we have calculated Polyakov-loop correlation 
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functions for 3-flavour QCD at a temperature close but below the transition tem- 
perature, T ~ 0.97 T c . This analysis has been performed with the p4-action on 
lattices of size 16 3 x 4 as a by-product of our determination of the critical couplings 
at various quark mass values. In addition we used the string tensions obtained from 
our zero temperature calculations at a nearby value of the gauge coupling, corre- 
sponding to T = T c , to set the scale for the potential and the spatial separation r of 
the static quark anti-quark pair. These potentials are shown in Figure [5|. They have 
been normalized such that they agree with the zero temperature Cornell potential 
at short distances, i.e. at r = 1/4T. 

We note that with decreasing quark mass the screening effects resulting from 
the presence of virtual quark anti-quark pairs lead to a more rapid flattening of the 
potential. While for m q = 1.0 the potential still coincides with the finite temperature 
potential in the pure gauge sector |36j at least up to r^/a = 3, the influence of 



dynamical quarks becomes large already for m q = 0.6. At small quark masses the 
screened potential reaches a limiting form and starts deviating from the pure gauge 
heavy quark potential already at r^fo ~ 0.7 and saturates at r^fo ~ 1.5. This 
corresponds to distances r ~ 0.3 fm and 0.7 fm, respectively. In general we find 
that deviations from the chiral limit are small for quark masses below m q ~ 0.1 or 
mps — 1.8-s/a. 

The above analysis of the quark mass dependence has been performed at a fixed 
temperature close to T c . In addition screening effects are, of course, also strongly 
temperature dependent. In the pure gauge case this leads to a temperature depen- 
dence of the string tension |3l| . In the presence of dynamical quarks we find that 
string breaking sets in at larger distances when the temperature is lowered. Note 
that this happens although at fixed bare quark mass m q the physical quark mass 
decreases with decreasing temperature. Thermal effects seem to become small for 
temperatures below T ~ 0.65 T c . This becomes evident from Figure [] where we 
show results from a calculation for 3-flavour QCD with quark mass m q = 0.1. In 
the low temperature regime the potential starts deviating from simple Cornell-type 
confinement potential for r^fo ~ 1.2 and reaches a plateau for r^fa ~ 3. This 
corresponds to distances r ~ 0.6 fm and 1.4 fm, respectively. 

To quantify the screening effect in the presence of light dynamical quarks we 
define the depth of the potential as the difference between the asymptotic value at 
large distances and the potential at distance r^fo = 0.5, 

AV = lim V(r) - V(Q.5/y/a) (4.3) 

which is the point at which the Cornell-potential vanishes if one chooses the cou- 
pling of the Coulomb term as a = 0.25. It approximately corresponds to a dis- 
tance r ~ 0.23 fm. In Figure |7] we show this difference of potential energies. 
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Figure 7: Temperature dependence of the depth of the heavy quark potential, 
defined in Eq. [4.3| , calculated in 3-flavour QCD with quarks of mass m q = 0.1. In 
order to convert to physical units we have used = 425 GeV as input. 

We note that AV^ changes most rapidly in the vicinity of T c . At T c we obtain 
AV ~ 0.5yfa ~ 200 MeV which, in fact, is compatible with the depth of a simple 
unscreened Coulomb potential, V{r)/^/~a= — a/r^/a with a = 0.25. 



5 Conclusions 

We have analyzed the quark mass and flavour dependence of the critical temperature 
of QCD using improved gauge and staggered fermion actions. We find that in a large 
interval of quark mass values, which corresponds to zero-temperature pseudo-scalar 
meson masses mpg between 360 MeV and about 2.3 GeV the transition temperature 
T c (mps) depends only weakly on mps- In this entire mass interval the transition 
temperature in two-flavour QCD is about 10% larger than that of three-flavour 
QCD. An extrapolation of both pseudo-critical temperatures to the chiral limit yields 
T c (0) = (173 ± 8) MeV in the case of two-flavour QCD and T c (0) = (154 ± 8) MeV 
for three-flavour QCD. In both cases we have quoted statistical errors only. From a 
comparison with calculations performed with unimproved actions and the experience 
gained with discretization errors in that well as in the pure gauge sector we 

estimate that systematic errors inherent in the present analysis are of the same order 
of magnitude. The results obtained here with staggered fermions are also consistent 
with results obtained from studies of two-flavour QCD with clover-improved Wilson 
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fermions 



Our analysis of the temperature and quark mass dependence of the heavy quark 
free energy shows that screening at T ~ 0.97T C is very efficient already in a quark 
mass regime where the (zero-temperature) pseudo-scalar meson masses are about 
1 GeV. Furthermore, screening effects become approximately quark mass indepen- 
dent for smaller masses. For a fixed value of the bare quark mass, m q = 0.1, we 
find that string breaking effects set in at a distance of about 0.6 fm for small tem- 
peratures, T < 0.6T C . With increasing temperature the onset point gradually shifts 
to shorter distances. At T c the heavy quark free energy gets screened already at 
distances of about 0.3 fm. 
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